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We introduce an extended version of oxDNA, a coarse-grained model of DNA designed to capture the ther¬ 
modynamic, structural and mechanical properties of single- and double-stranded DNA. By including explicit 
major and minor grooves, and by slightly modifying the coaxial stacking and backbone-backbone interactions, 
we improve the ability of the model to treat large (kilobase-pair) structures such as DNA origami which are 
sensitive to these geometric features. Further, we extend the model, which was previously parameterised to 
just one salt concentration ([Na~'“] = 0.5 M), so that it can be used for a range of salt concentrations including 
those corresponding to physiological conditions. Finally, we use new experimental data to parameterise the 
oxDNA potential so that consecutive adenine bases stack with a different strength to consecutive thymine 
bases, a feature which allows a more accurate treatment of systems where the flexibility of single-stranded 
regions is important. We illustrate the new possibilities opened up by the updated model, oxDNA2, by pre¬ 
senting results from simulations of the structure of large DNA objects and by using the model to investigate 
some salt-dependent properties of DNA. 


I. INTRODUCTION 

Deoxyribonucleic acid (DNA) performs the crucial 
function of storing genetic information in living organ¬ 
isms. It is made up of repeating units called nucleotides, 
each of which consists of a sugar and phosphate backbone 
plus a base (either Adenine (A), Guanine (G), Thymine 
(T), or Gytosine (G)) attached to the sugar. Watson- 
Crick base-pairing, A with T and G with G, along with 
planar stacking interactions between bases and the con¬ 
straints of the backbone, leads to the formation of the 
well-known double-helical structure of DNA. 

The Watson-Crick complementarity of DNA also per¬ 
mits the rational design of DNA objects for which the 
intended structure is the global free-energy minimum, a 
property which has been exploited to create a wide vari¬ 
ety of 2D and 3D nanostructures.^ilt^ Further, these DNA 
objects can be functionalised,®^ with potential applica¬ 
tions ranging from nanomedicine to nanoelectronics .1^ 

Theoretical and computational approaches to mod¬ 
elling DNA have been widely exploited to probe the 
behaviour of the molecule in both a biological and a 
nanotechnological context. At the finest level of de¬ 
tail, quantum chemistry calculations have been used to 
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study the interactions between nucleotidesalthough 
the high computational cost of such an approach limits 
these methods to interactions between nearest-neighbour 
base pairs in vacuum. Glassical all-atom approaches, 
where every atom of DNA and the surrounding solvent 
is modelled as a point particle with effective interac¬ 
tions, have been widely employed to study small DNA 
motifs,®®land have recently been applied to larger DNA 
systems.^i®^^ However, simulating rare-event processes 
such as the breaking or formation of base pairs remains 
a challenge with these models, with ps time scales being 
the limit of what is currently accessible.®^ At the other 
end of the scale, theoretical approaches have been de¬ 
veloped to understand certain large-scale properties of 
DNA. These include the wormlike-chain mod^ which 
treats DNA as a continuously flexible polymer.®! While 
such models can provide useful insights into the physi¬ 
cal properties of DNA, they are not detailed enough, by 
design, to address processes such as duplex formation. 

The middle ground between analytical and all-atom 
approaches is occupied by coarse-grained DNA models. 
Such models integrate out many of the degrees of free¬ 
dom of the DNA nucleotide, and often neglect the solvent 
molecules; these approximations inevitably imply a com¬ 
promise between accuracy, generality and computational 
efficiency, so that care must be taken in applying these 
models to a given problem. However, the simplified pic¬ 
ture presented by such an approach can be a strength, as. 
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in addition to greatly increasing the time scale and num¬ 
ber of nucleotides that can be studied, it can allow one 
to understand the generic physics governing the system 
more easily. 

Many coarse-grained DNA models have been devel¬ 
oped in recent years^^^^ in this work w e im prove 
the model due to Ouldridge et ai, oxDNA.l23^ This 
nucleotide-level model is designed with a heuristic, “top- 
down” approach, with a focus on reproducing well-known 
properties of DNA (such as the helical structure of the 
B-DNA duplex) and experimental results (such as du¬ 
plex melting temperatures), rather than, for example, 
building the model up by integrating out details from 
an all-atom representation. In the oxDNA model, each 
nucleotide is a rigid body with three interaction sites 
that have mutual, highly anisotropic interactions. This 
treatment is sufficient to obtain good agreement with ex¬ 
perimental data on the structural, mechanical and es¬ 
pecially the thermodynamic properties of single- and 
double-stranded DNA. Consequently, the model has pro¬ 
vided key insights into many different processes rele¬ 
vant to DNA nanotechnolog}i2lHlZl and biophysics ,l2SlE2l 

and importantly, has also been shown to provide direct 
agreement with experimentally measured properties on 
a range of systems including DNA overstretching,!^ a 
two-footed D NA w alker,!^ and toehold-mediated strand 
displacement EHm 

Despite these achievements, there are some areas 
where oxDNA can be improved. The model was pa- 
rameterised to [Na’*’] = 0.5 M, a high salt concentra¬ 
tion similar to that used for many applications in DNA 
nanotechnology. However, the ability to study DNA be¬ 
haviour as a function of salt would allow quantitative 
comparison with a greatly expanded set of experiments, 
and, in particular, if we wish to apply oxDNA to biolog¬ 
ical systems we would like to work at physiological salt 
([Na'*’] « 0.15M). At the same time, the wealth of ex¬ 
perimental data for the thermodynamic and mechanical 
properties of DNA as a function of salt concentratioiP^^!^ 
makes fitting the salt-dependent properties of an ex¬ 
tended model possible. 

While the detailed e ffect o f salt on DNA electrostatics 
can be highly complex,l2SllIl here we use a simple Debye- 
Hiickel interaction term, as first implemented for oxDNA 
by Wang and Pettitt in Ref. [39l to model how salt screens 
repulsive interactions. The more coarse-grained descrip¬ 
tion is commensurate with the level of approximation 
used more generally in the oxDNA framework. We care¬ 
fully parameterise the model to an extended set of exper¬ 
imental data for melting temperatures and persistence 
lengths, including the behaviour of single strands. 

A second area that merits attention is the perfor¬ 
mance of the model in simulating the structure of large 
(kilobase-pair) DNA objects. The model value of the 
B-DNA pitch did not come under too much scrutiny 
in the original parameterisation as there is some dis¬ 
agreement in th e liter ature about the precise value of 
the average pitch.!^^!^®^ However, improved experimen¬ 


tal tec hnique s in fabrication and imaging of large DNA 
objectjn^SES] have presented an opportunity to finely 
tune this value, because small adjustments to the duplex 
pitch can result in signihcant changes to the global twist 
of a large-scale DNA nanostructure. In addition, simu¬ 
lating these large-scale structures has illustrated the po¬ 
tential importance of nicks and junctions for the effective 
duplex pitch, so that improving the model’s description 
of these effects has become a priority. Finally, the origi¬ 
nal model duplex had grooves with equal widths, whereas 
B-DNA is known to have a larger major groove and a 
smaller minor groove. This implies that the positions of 
the backbone sites in the model, which are directly re¬ 
lated to the groove widths, could be more realistic. This 
detail could be relevant, for example, in origami struc¬ 
ture, where the precise backbone positions are known to 
be important for junction placement.!^ 

The oxDNA model was previously given sequence- 
dependent hydrogen-bonding and stacking strength^^by 
fitting to the duplex thermodynamics of the SantaLu- 
cia model.!^ As the SantaLucia model gives results at 
the base-pair step level, one can only extract the aver¬ 
age stacking strength for the two stacking interactions 
present in a given base-pair step. In particular, this 
means that in oxDNA the AA and TT stacking interac¬ 
tions are the same, whereas it is well known that the AA 
stacking interaction is significantly stronger than the TT 
interaction, an important property for DNA nanotechnol¬ 
ogy, for example, where poly-T sing le-stranded regions 
are often used as flexible linkers To remedy this we 
use new experimental data to reparameterise the AA and 
TT stacking interactions in the model. 

In the following sections, after briefly introducing the 
original oxDNA model, we consider each change to the 
model in turn, and then we highlight some important 
aspects of the behaviour of the new model. 

II. THE ORIGINAL OXDNA MODEL 

In the oxDNA model introduced by Ouldridge et 
strands of DNA are represented by a chain of rigid bod¬ 
ies, with each rigid body representing a nucleotide. The 
coarse-grained, pairwise potential for the model can be 
written as a sum of interaction terms 

V = (14>ackbone + 14tack + ^exc) + 

nearest neighbours 

(VhB + 141.088 stack T I4xc T I4oax stack)- 

other pairs 

The nearest neighbours (adjacent nucleotides on a 
DNA strand) interact via ^backbone, 14tack, and t44c: 
which represent the connectivity between neighbouring 
backbones, the favourable stacking interactions between 
neighbouring bases, and excluded volume terms, respec¬ 
tively. All other nucleotide pairs interact with Vhb, 
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FIG. 1. Schematics contrasting the original oxDNA model 
(left), with equal groove widths, with oxDNA2 (right), which 
has differentiated major and minor grooves, (a) A cross sec¬ 
tion of a duplex with one base pair displayed. The large 
dashed circle shows the helix radius, and each nucleotide is 
represented by three circles joined by a line; the large solid cir¬ 
cles represent the backbone sites, while the small solid circles 
represent the stacking (closer to the backbone) and hydrogen¬ 
bonding (at the end of the nucleotide) sites. For oxDNA2, a 
value of 20° was chosen for the angle 7 . (b) A representation 
of a DNA duplex for each model. 


^ross stack: ^xc ctucl Vcoax stack: Corresponding to hy 
drogen bonding between complementary bases, cross¬ 
stacking, excluded volume, and coaxial stacking between 
non-nearest neighbours, respectively. We now highlight 
aspects of the original model that are relevant to the im¬ 
provements made in this paper. 

Firstly, in its original formulation, the model was pa- 
rameterised for a sodium ion concentration of 0.5 M, 
chosen to reflect the high salt concentrations commonly 
used in DNA nanotechnology. The electrostatic interac¬ 
tions between the negatively charged phosphate groups 
on the DNA backbone were incorporated into the back¬ 
bone site’s excluded volume, an approximation which can 
be justified by the very short Debye screening length at 
that relatively high ion concentration. 

Secondly, the original oxDNA modeP^ represented 
each nucleotide as a linear rigid body (Fig. [^. The 
optimal configuration for base-pairing occurs when the 
two nucleotides point directly at each other. As a conse¬ 
quence, the DNA double helix was symmetric, with the 
two grooves having equal widths. 

Thirdly, in the original oxDNA model introduced in 
Ref. Ea all four types of base were treated equally ex¬ 
cept that only A-T and G-C base pairs could be formed. 
Later Sulc et alWi introduced sequence-dependent ther¬ 


modynamics into the model by making the strengths 
of the hydrogen-bonding and stacking terms depend on 
the identities of the interacting nucleotides. The nearest 
neighbour DNA model of SantaLucia,l^to which oxDNA 
was parameterised, does not resolve the difference be¬ 
tween AA and TT stacking as it works on the base-pair 
step level - therefore AA and TT stacking strengths were 
set to be the same in Ref. [571 

In this work we mostly work from the original, 
sequence-averaged parameterisation of the model rather 
than the sequence-dependent one, as it is more effi¬ 
cient to fit the thermodynamic parameters to sequence- 
averaged duplex melting temperatures as given by the 
SantaLucia model. The exception is the parameterisa¬ 
tion of the AA and TT stacking strengths, which did 
use the sequence-dependent parameters from Ref. |57| as 
a starting point, to allow the best possible comparison 
between the model and the experimental results that 
were used for the htting. After the parameters for 
oxDNA2 had been obtained, including new values for 
the sequence-averaged hydrogen bonding and stacking 
strengths, we then rescaled the sequence-dependent in¬ 
teraction strengths from Ref. |57| accordingly for use with 
the new model. 


III. INTRODUCING DIFFERENT WIDTHS FOR MAJOR 
AND MINOR DNA GROOVES 

B-DNA in the original oxDNA model has equal groove 
widths, while in reality DNA has a larger major groove 
and a smaller minor groove. Having realistic widths for 
the major and minor grooves is equivalent to having ap¬ 
propriately positioned backbone sites in the model, an 
important feature for the physical properties of many 
DNA motifs. For example, in DNA origami antiparal¬ 
lel double helices are joined by crossovers, for which the 
position of the ba ckbo ne has been shown to be crucial for 
origami structure.^^^ Another example is anisotropic du¬ 
plex bending: the duplex can be expected to bend more 
easily into the major groove than into the minor, if the 
groove widths are unequal. 

The oxDNA nucleotide is composed of three interac¬ 
tion sites: the hydrogen-bonding, stacking and backbone 
sites. We introduce different groove widths by changing 
the position of the backbone site while keeping the duplex 
radius unchanged (Fig.[^, such that, rather than lying on 
a straight line, the three interaction sites lie in a plane. 
The new nucleotide shape introduces an additional pa¬ 
rameter into the model, the angle 7 between the line 
from the duplex centre to the backbone site and the line 
from the duplex centre to the stacking site (Fig. Ba))- 
Given the coarse graining of the 18 atoms of the sugar- 
phosphate DNA backbone into a single interaction site, 
there is no definitive choice for the precise position of the 
backbone site and thus the value of the model parame¬ 
ter 7 (Fig. Ba)). We set 7 = 20°, a value which maps 
onto a full-atom representation of a DNA duplex well by 
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visual inspection, although values of 7 between 15° and 
25° would give an equally satisfying visual match. 

The backbone site is moved such that the duplex ra¬ 
dius is unchanged, and we note that the modification has 
a negligible added computational cost when simulating 
the model. However, the thermodynamic and mechan¬ 
ical properties are slightly affected. For the thermody¬ 
namics, we found a change of 1-2 K in the duplex melt¬ 
ing temperatures, and we modified the hydrogen-bonding 
and stacking strengths using the histogram reweighting 
method described in Section II C 2 of Ref. EH so that 
the agreement with experimental melting temperatures 
was as good as for the original model. The mechanical 
properties of DNA are less well constrained experimen¬ 
tally and so were not refitted. The mechanical properties 
for the new model can be found in Section IVlIIl 

One illustration of the importance of the groove widths 
in oxDNA for the structural properties of DNA assem¬ 
blies is provided by systems of 3-arm star tiles that are 
designed to form triangular prismatic polyhedra. We find 
that modifying the groove widths qualitatively changes 
the structure of trimers of these tiles (Fig. |^. Specif¬ 
ically, the body of the trimer defines a plane with two 
distinct faces. Zhang et al^ found that one of two possi¬ 
ble isomers of the polyhedron preferentially formed, im¬ 
plying that the free arms of the trimer systematically 
pointed in the direction of one of these two faces. We 
find a consistent result when the groove widths specified 
by oxDNA2 are used. When equal-sized grooves are used 
(as in the original oxDNA model) the trimer arms point 
in the opposite direction. 


IV. EFFECTIVE ELECTROSTATIC INTERACTIONS 

One major improvement presented in this paper is the 
introduction of a salt-dependent interaction term in the 
oxDNA model. Since the original goal of the oxDNA 
model was to simulate nanotechnology experiments, the 
thermodynamic and structural parameterisation was car¬ 
ried out at high salt. The very short electrostatic screen¬ 
ing length at these conditions allows one to incorporate 
the electrostatics into a soft excluded volume, somewhat 
circumventing the necessity of a proper treatment. The 
original parameterisation of the thermodynamics was 
carried out at 0.5 M [Na+J, a high enough value that 
further increasing it does not significantly change the 
physics at our level of coarse-graining. 

The problem of treating electrostatics properly for 
DNA in solution is a complicated one. Perhaps the most 
evident issue is that the typical dimensions of nucleotides 
are comparable to the Debye length of the solution, ren¬ 
dering a mean-field treatment hard to justify. Also, in 
some cases the presence of salt ions affects the local 
structure of nucleic acids, by stabilizing some arrange¬ 
ments or destabilizing others. Ion condensation may also 
lead to stronger screening of the electrostatic interactions 
than otherwise expected - this has been incorporated 


into co arse-g rained DNA models through partial effective 
chargesThus, in principle, many non-trivial effects 
must be taken into account when modelling the electro¬ 
static interactions, and the de bate o n the best way to do 
so implicitly is still unresolved.l^^^^ 

Here, we choose a very simple treatment, based on the 
Debye-Hiickel model for screened electrostatics. This ap¬ 
proach has b een used previously for other coarse-grained 
models,l23|63| and was first introduced into oxDNA by 
Wang and Pettitt.^^We note that this treatment is con¬ 
sistent with the coarse-grained nature of the model, and 
we use the same top-down strategy that was used in 
the original parameterisation to design the effective elec¬ 
trostatic interactions: our goal is to introduce a term 
in the potential that will reproduce the thermodynamic 
and mechanical effects of salt concentration on DNA, and 
thus should be regarded as an effective interaction rather 
than an attempt to rigorously model the local effects of 
charges. This should be kept in mind when interpreting 
the results obtained with oxDNA, in particular at low 
salt concentration. 

Since the modelling of electrostatics is rather crude, 
we restrict our parameterisation to salt concentrations of 
O.IM of monovalent salt or greater. This restriction is 
also due to the fact that we parameterise our thermody¬ 
namics to the model of SantaLuda,!^ which was fitted in 
a similar salt regime. Importantly, physiological condi¬ 
tions fall within this range, which will allow quantitative 
comparison between simulations of DNA systems with 
our model and experiments at physiological conditions. 

The detailed Debye-Hiickel form which is added to the 
non-bonded interactions in the potential of Eq. takes 
the following form: 

( 2 ) 

ij U 

where g'eff is the effective charge situated at the backbone 
site of each of the nucleotides, is the magnitude of 
the distance between the backbone sites of nucleotides 
i and j, cq is the permittivity of the vacuum, is the 
relative permittivity of water, e is the elementary charge 
and ( 7 eff is a dimensionless effective charge. In principle, 
tr depends on Vij , 1 ^ and weakly depends on temperature 
and salt concentration. However, for oxDNA2 we set 
Cr to be a constant value, in keeping with the coarse¬ 
grained approach taken for the rest of the model. In 
particular, we choose = 80, the standard value for 
water. In Eq. we have stressed that the interaction 
depends on the temperature T and on the (monovalent) 
salt concentration I through the Debye length Adh(T, I): 


Adh(T,/) = 


eoCr/CBT 

2Np^e^I ’ 


(3) 


where Aa is Avogadro’s number and fee is Boltzmann’s 
constant. 

To ensure the computational efficiency of simulating 
the model, we set the interaction to zero at a finite dis¬ 
tance; and to allow simulation with molecular dynamics. 
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FIG. 2. In oxDNA, the structure of the DNA trimer due to Zhang et alW^ is found to qualitatively change as a function of 
duplex groove angle 7 . A representative trimer configuration with 7 = 20° is shown in (a) from the side and with 9, the angle 
between each arm and the plane of the main, triangular trimer section, displayed and in (b) from the top. (c) Shows the average 
value of 0 as a function of 7 , with 7 = 0° corresponding to equal groove widths as for the original oxDNA model, and 7 = 20° 
being chosen for oxDNA2. 


we introduce a quadratic smoothing potential so that the 
interaction goes to zero smoothly. The quadratic smooth¬ 
ing, the details of which are reported in Appendix |A 1[ 
is introduced after a cutoff rsmooth, which we choose to 
be 3 Adh- This cutoff allows us to use all the standard 
techniques to improve the simulation efficiency via the 
use of Verlet lists and/or cells.We have checked that 
introducing our chosen cutoff, rsmooth = 3 Adh , has a neg¬ 
ligible effect on the duplex thermodynamics results used 
to parameterise the interaction (Fig. [^a)). 

Our representation of DNA uses a single rigid body 
per nucleotide, and the best choice of where to put the 
charge is not obvious. All the atoms of the sugar and 
phosphate groups of the backbone are represented in a 
single interaction site, and it is thus natural to put the 
charge, which in real DNA is located on the phosphate, 
on that interaction site. Importantly, the backbone site 
of a nucleotide is placed almost in between the phosphate 
of that nucleotide and the phosphate of the neighbouring 
one, which could potentially lead to some unphysical ef¬ 
fects. Also, we should stress that having a charge at each 
backbone site means that the DNA has as many charges 
as nucleotides, which is not always true in real systems: 
very often, the terminal phosphate at the 3' end is cut off, 
removing a charge. The absence of this charge can cause 
measurable effects on the thermodynamics, and indeed 
the SantaLucia modeP^ requires the presence or absence 
of the terminal phosphates as an input parameter. In 
keeping with our coarse-graining approach, we put a half 
effective charge on each of the terminal nucleotides to in¬ 
corporate the fact that each charge should be halfway in 
between our backbone sites, emulating a system with the 
terminal charge removed, and parameterise to the San¬ 
taLucia model in a way consistent with this approach. 

The parameter that we have tuned to reproduce the 
thermodynamics predicted by the model of SantaLucia 
is the effective charge geS- To do this, we used thermo¬ 
dynamic integration^ to compute the melting tempera¬ 
tures of duplexes of length 5, 6, 7, 8, 10, 12 and 15 as a 


function of Qes at several salt concentrations, and chose 
the value that best reproduced the melting temperatures 
predicted by SantaLucia’s model. 

The melting temperature of a duplex is defined as the 
temperature at which half of the strands in a stoichio¬ 
metric bulk solution are in the duplex state. We cannot 
simulate the bulk system with oxDNA. However, at the 
bulk melting temperature for a given concentration of 
strands, the bound and unbound free energies obtained 
for a system of two strands at the same concentration, 
Fb and Fub, are related byP^ 

FhiTm, qes) = T'ub(T'm, Qes) — ln(2), (4) 

where we have made explicit the dependence of F on 
the effective charge geff- The constant on the right hand 
side of Eq. accounts for concentration fluctuations that 
are present in bulk but not for two strands in a periodic 
box.ES The equation is exact in the dilute limit, which is 
indeed where the thermodynamics of duplex formation 
are usually studied and is the relevant limit for oxDNA. 

We can take advantage of the relation of the free en¬ 
ergies given by Eq. [^for the thermodynamic integration. 
Eor small changes in T and ges, it is possible to write 


F(Tm -I- dF, geff -I- dgeff) = 
dF dF 

F{T^, geff) -H ^dT -h - -dgeff = 

oT dgeff 


geff) -I- 


F-{V) 


dV 

W 


dT + 2^^dgefr. 


Qeff 


(5) 


The last step includes the unusual term {dV/dT). This is 
because our potential, like many coarse-grained models, 
depends explicitly on T through 14tack and Vdh- 

To compute the change in Fn introduced by a small 
change in geff, one can impose the condition in Eq. [^at 
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FIG. 3. (a) Melting temperature of duplex DNA of differ¬ 
ent lengths I at [Na''']=0.2M as a function of the effective 
charge q^a obtained with thermodynamic integration, (b) 
Average melting temperature difference between our model 
and the SantaLucia modeP^ for different values of 1 as a 
function of geff. Each point corresponds to the average over 
[Na''']=0.1, 0.2, 0.3, 0.4, 0.5M of the magnitude of the dif¬ 
ference in melting temperature between our model and the 
SantaLucia model. The plot suggests an optimal value for 
qeff of 0.815, indicated by the vertical dashed line. 

the new 

Fh{Tyn + dr, qes + dgeff) = 

Fuh{T^ + dr, qeff + dqeff) - k^{T^ + dr) ln(2), (6) 

and, by using Eqs. [^and|^ one obtains 
drm _ (2/geff)((VDH)ub — (VDH)b) 

^ " ((Kb) - {Vf,))/T^ - ((9Eub/9r) - {dVf,/dT)y 

(7) 

Eq. is a differential equation that allows us to follow 
the change in melting temperature as geff is changed. 
We note that Eq. is an extended Clausius-Clapeyron 
relation,!^ and the quantities on the right-hand side are 
readily accessible with separate simulations of the bound 
and unbound states. Since we are dealing with small sys¬ 
tems, all these simulations are very quick and it is thus 
easy to achieve very accurate results. As a starting point 


for the thermodynamic integration we use the melting 
temperatures from SantaLucia, which our original model 
(equivalent to the current model at q^ff = 0) reproduces 
within tenths of a Kelvin. Since both Ktack and Vdh only 
weakly depend on temperature, we used Eq. [^assuming 
{dV^h/dT) — {dVh/dT) = 0 to obtain an optimal value 
for qeff, and checked a posteriori the results of thermo¬ 
dynamic integration with melting simulations. We have 
found that the melting temperature differences computed 
with Eq. [^are accurate to within 0.2 K. 

For a range of salt concentrations and duplex lengths, 
we performed thermodynamic integration to find the 
melting temperature as a function of qeff. Some of the 
results are depicted in Fig. The melting temperature 
of the duplex decreases with increasing effective charge 
qeff, but not dramatically so (Fig. [^a)). This is be¬ 
cause the increased interchain repulsion in the duplex 
state is partially balanced by the lower entropy of the 
single-stranded state (due to it adopting a slightly more 
extended state to reduce intrachain repulsion). The data 
in Fig. [^b) show that the qeff at which the difference in 
melting temperature between oxDNA and SantaLucia is 
minimised depends on the duplex length and lies in the 
range 0.75 < qeff < 0.95. The best overall predictions are 
obtained for an effective charge qeff = 0.815. It is reassur¬ 
ing that this value does not deviate significantly from 1 
(corresponding to the value given by Debye-Hiickel the¬ 
ory), and that the best value for qeff varies little with 
duplex length. In addition, it is not uncommon t o use 
a value of qeff < 1 for coarse-grained DNA models.^S^Illl 
One argument is that ion condensation, which is known 
to occur for DNA, will screen the phosphate charges more 
strongly than expected from Debye-Hiickel theory. This 
will lead to a lower effective charge when fitting a model 
using a Debye-Hiickel treatment to experimental results, 
although such arguments should be applied with caution 
to a crude mean-field approach such as this one. 

We note that introducing this explicit electrostatic 
term in our model potential will raise the computational 
expense of each simulation step compared to the original 
oxDNA model, as the electrostatic term will generally re¬ 
sult in oxDNA2 having a larger interaction range than the 
original model has. This effect increases at lower salt, as 
the Debye-Hiickel term becomes more long-ranged. For 
example, we find that simulating a 10-bp duplex with 
MD for a given number of steps takes 1.4x as long with 
oxDNA2 at [Na''']=0.5M and 1.7x as long with oxDNA2 
at [Na“'']=0.1 M as it did for the original oxDNA. 


V. IMPROVING STRUCTURE PREDICTION FOR 
LARGE-SCALE DNA OBJECTS 

The model was originally parameterised for small 
single- and double-stranded DNA structures. If we wish 
to study large-scale structures, we need to ensure that we 
can reproduce the existing experimental data for these 
larger constructs. A good test case is provided by the 
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FIG. 4. Simulation snapshots for the L- (left), N- (middle), 
and R- (right) type helix bundles, each composed of roughly 
15,000 nucleotides. The designs are taken from Dietz et al.l^In 
experiment and in oxDNA2, the N-type helix bundle exhibits 
close to zero global twist. 


Global Twist/° 

Helix Bundle Type oxDNA2 experiment 


N 

0.2 

0 

L 

-22.2 

-31 +/- 5 

R 

26.9 

26 +/- 5 


TABLE I. Global twist in the DNA helix bundles of Dietz 
et for simulations using oxDNA2 and experiment. The 
oxDNA2 simulations were run with [Na+]=0.5 M for 1.8 x 10® 
molecular dynamics (MD) steps for each design. 


work of Dietz et aZ.,li^who measured the global twist of 
three different DNA origami structures, described as 10- 
by-6 helix bundles (Fig. |^. We denote the three origami 
structures as L-, N- and R-type. In the experiment, 
they were designed to impose a pitch of 10, 10.5, and 11 
base pairs per turn (bp/turn) on the constituent DNA 
double helices and these different designs exhibited left- 
handed, no, and right-handed global twist, respectively, 
when multimerised to form ribbons and visualised with 
transmission electron microscopy. One might think that 
this result implies that DNA has a natural pitch of 10.5 
bp/turn, since the design with that inherent periodicity 
did not result in a globally twisted system. However, it 
is not that straightforward. 

Unsurprisingly, when simulated with the original ver¬ 
sion of oxDNA, the N-type origami still showed a signifi¬ 
cant right-handed global twist, chiefly because the duplex 
pitch for the model was 10.36 bp/turn. However, even 
when simulated with a version of oxDNA modified so 
that the model duplex pitch was 10.5 bp/turn, the N-type 
origami displayed a right-handed twist. One reason for 
this is the following: in oxDNA, the helical twist across 
nicks and junctions is larger by about 3.9° and 2.5° de¬ 
grees respectively than for a normal duplex step (we call 
this an overtwist), so that fractionally fewer base pairs 
are required for a helix turn than would otherwise be ex¬ 
pected. Although these differences are small, they add 
up constructively to create a global twist on a structure 
that would otherwise have no such twist. For a nicked 
duplex in oxDNA, the overtwist occurs because, oppo¬ 


site the nick, the non-nicked strand prefers a larger twist 
than for duplex DNA (just as stacked single strands do 
in the model), and the nicked backbones lack the FENE 
spring that would usually oppose this tendency. A simi¬ 
lar argument applies for strands with junctions. 

Although there is no direct evidence available to show 
whether this overtwist is physically realistic, there are 
multiple lines of experimental evidence that suggest that 
the pitch of duplex DNA is close to 10.5 bp/turn.E3SSl 
Given the evidence from Dietz et that the effective 
pitch in DNA origami structures is also close to 10.5 
bp/turn, we decided to reduce the overtwisting in oxDNA 
as much as is possible. To achieve this, we modify the 
coaxial stacking term of the potential, Vcoax stack, so that 
the overt wist i s 0° for a junction and 1.3° for a nick (see 
Appendix | A 3| for details of the changes to the potential). 
We have verified that this change has very little effect on 
the other features of the model. 

Even when the overtwist at nicks and junctions has 
largely been removed, we require an intrinsic duplex pitch 
of 10.55 bp/turn in order for the oxDNA N-type helix 
bundle to have zero global twist, and we set the pitch 
to this value (at [Na“'']=0.5 M) for oxDNA2 by modify¬ 
ing Vbackbone- We found that this modification to the 
backbone potential changed the duplex melting temper¬ 
atures by around 2 K, and we refitted the model’s hydro¬ 
gen bonding and stacking strengths to correct for this us¬ 
ing the same method as was used for the thermodynam¬ 
ics refitting when implementing unequal helical groove 
widths. The requirement for a pitch of 10.55 bp/turn, 
rather than the 10.5 bp/turn that one might expect, is 
due to a few subtle effects that are currently being inves¬ 
tigated further. 

The global twist measured for the helix bundles with 
the new oxDNA2 model (Appendix B 1| describes how 


the global twist was measured from simulations) is com¬ 
pared to the experimental results in Table B while typi¬ 
cal simulation snapshots are shown in Fig. ^ The slight 
modifications to Vbackbone used to set the model pitch are 
given in Appendix |A 2| We note that the modifications 
to Vbackbone changed the duplex melting temperatures in 
the model by 1-2 K, which were then refitted using his¬ 
togram reweighting (as described in Section H C 2 of 
Ref. I61|l to give an agreement with experimental melting 
temperatures that was as good as for the original model. 


VI. AA/TT SEQUENCE DEPENDENCE 

As a further improvement, the oxDNA2 model in¬ 
corporates a more realistic sequence-dependent stack¬ 
ing interaction, which is achieved by differentiating be¬ 
tween the AA and TT stacking interaction strengths. 
In the previous parameterisation of oxDNA,the 
sequence-dependent base-pairing and stacking interac¬ 
tion strengths were obtained by fitting the oxDNA du¬ 
plex melting temperatures to the SantaLucia model, a 
nearest-neighbor model that is able to predict experimen- 













tal duplex melting temperatures very well.^ The San- 
taLucia model is designed at the level of base-pair steps, 
where each base-pair step consists of four bases, with a 
free-energy difference between its single-stranded state 
and its duplex state. In total there are 10 unique base- 
pair steps in the SantaLucia model: AA/TT, AT/AT, 
TA/TA, GC/GC, GG/GG, GG/GG, GA/TG, AG/GT, 
TG/GA, GT/AC (for example, AG/GT refers to a base- 
pair step with complementary bases AG on one strand 
and CT on the other, both specified in the 3' to 5' di¬ 
rection). Therefore, fitting oxDNA to the SantaLucia 
model only allows one to find the sum of the strengths 
of the two stacking interactions between the nucleotides 
that are within a base-pair step. For example, it is only 
possible to find the average of the AA and TT stacking 
strengths from the SantaLucia model. 

However, experimental evidence from sequence- 
dependent measurem ents o f the mechanical properties of 
single-stranded DNA,I^ ^^^ as well as hairpin stabilities 
and closing rates ,122122 have revealed that sequences of A 
bases are much stiffer than equal length sequences of T 
bases. It has been argued that this is evidence that con¬ 
secutive AA bases stack much more strongly than TT 
bases do. Here we use original experimental data on the 
stabilities of hairpins with either a poly-A or a poly-T 
loop, to differentiate between the AA and TT stacking 
strengths. All hairpins have 6-bp stems and loops of ei¬ 
ther 21 or 31 bases. Details of the experimental systems 
and the sequences are given in Appendix [C| 

To parameterise our model to reproduce the experi¬ 
mental data, our procedure is to vary the AA stacking 
interaction strength, caa, while fixing the sum of the AA 
and TT stacking interaction strengths to the value that 
is obtained by fitting oxDNA’s duplex melting tempera¬ 
tures to the SantaLucia model, i.e. 


EAA + CTT — 2eavgi 


( 8 ) 


tegrals 


F^{I,x,qes) = Fi°'> + 


, 8 ) 


Si 


eff 


dx' 


fdFfyi^I^X ,Qeff) 
dx' 


where x measures the deviation of the AA and TT 
strengths from Cavg, such that Eq. [^is satisfied, i.e. 

caa = (1 + a:)eavg, ctt = (1 - a:)eavg • (10) 

Taking the derivatives of the free energy in Eq. we 
obtain, 


F, 


Jo JeS 


dx' 


1 -I- X' 


f (kstack( 2^5 ^ ; Qeff)) 


AA 


1 — X' 


f (Lstack(d, X , Qeff )) 


TT 


( 11 ) 


where ( 14 tack( 1^5 x, is the average stacking en¬ 

ergy of all AA(TT) nucleotides in a strand in state a with 
salt concentration I, x dehned as in Eq. and back¬ 
bone charge q^s. The terms appearing in the integrands 
of Eq.[^can be obtained by running short simulations of 
bound and unbound states separately and therefore the 
free energies can be calculated in an efficient manner. 
The results for the difference in 5/S.F between oxDNA 
and experiment as a function of caa for the hairpins 
with 21- and 31-base loops are shown in Eig.[^ the value 
for the AA stacking strength that minimises this differ¬ 
ence is found to be caa ~ l.OTScavg (corresponding to 
ett ~ 0.925eavg), which is not too dissimilar to the pre¬ 
liminary value suggested in Ref. We note that this 
value gives satisfactory predictions for 5/S.F for a wide 
range of salt concentrations down to 0.05 M. 


where ctt is the TT stacking strength, Cavg is the 
strength for both AA and TT obtained from fit¬ 
ting oxDNA to the SantaLucia model. We vary 
the AA stacking interaction strength to match the 
experimental differences of the thermal stabilities of 
the poly-A-loop and poly-T-loop hairpins, specifically 
S^F — AF)a— loop) loop); where AF/x—loop) ~ 

■Fb,(x-ioop) - A/b,(x-ioop) is the difference in free en¬ 
ergy between the bound (hairpin) state and the unbound 
(open) state for a hairpin with a poly-X loop. 

We obtain the bound and unbound free energies Fa, 
where a € {b,ub}, by using thermodynamic integration. 
We start by calculating Fa\ the free energy of a refer¬ 
ence state for which the effective charge on the backbone 
site, 9eff, is set to 0, and eaa = ctt = Cavg- Fa{I, x, geff), 
the free energy of the state with salt concentration I, 
stacking strengths eaa and ctt; and effective backbone 
charge qes, is then obtained by solving the following in- 


VII. THE OXDNA2 MODEL 

In summary, the oxDNA2 potential can be written as 

KxDNA 2 = ^ (K*ackbone + Ktack + Kxc) + 

nearest neighbours 

(biB + b:ross stack + f4xc + i^oax stack + bon); 

other pairs 

( 12 ) 

where a V* indicates that the term is either modified 
for oxDNA2, or, in the case of VSh; aew in oxDNA2. 
The modified parameters for oxDNA2 are compared with 
those for oxDNA in Table |Al[ and a full account of the 
changes is given in Appendix All other parameters 
remain the same as in the original model. 













9 





FIG. 5. Deviation of oxDNA’s value for 5AF from experiment 
at 295.6K, as a function of CAA/tavg- The grey region contains 
the zero-deviation points for every curve. The experimental 
results are shown in Fig. |A5| and details of the experimental 
setup can be found in Appendix [Uj 


FIG. 6. (a) Rise and (b) pitch of a 60-bp duplex as a function 
of salt concentration in oxDNA2. For each salt concentration 
the duplex was simulated for at least 3 x 10® MD steps. The 
error bars, which are narrower than the plot markers, show 
the standard error on the mean given by averaging over 10 
independent estimates for each data point. 


We emphasise that, after all the relevant changes to 
the potential were made, the hydrogen bonding and 
stacking parameters were modihed to ensure that the 
close agreement to experimental duplex melting temper¬ 
atures achieved for the original model was retained with 
oxDNA2 (this was done immediately after the changes to 
^backbone, as described at the end of Section [v|). 


VIII. PHYSICAL PROPERTIES OF OXDNA2 

The structural, mechanical and thermodynamic prop¬ 
erties of DNA for the new version of oxDNA presented 
in this paper, which we call oxDNA2, are slightly differ¬ 
ent from the properties for the original oxDNA model. 
We briefly highlight the most important of these changes 
here; in addition these properties are given as a func¬ 
tion of salt concentration as this is now possible with the 
new model. We note that all of the results described 
in this section were computed using the final version of 
the oxDNA2 model as summarized in Section I VIII The 
details of the simulations used to compute these results 
are given briefly in the figure captions and in detail in 
Appendix |B 3| and Table [X3l 

The structural properties of the new model, specif¬ 
ically the pitch and rise of double-stranded DNA, are 
presented in Fig. As might be expected, the rise in¬ 
creases with decreasing salt concentration, due to the 
greater repulsion between backbone sites. The pitch also 
increases with decreasing salt, consistent with the mea¬ 
sured increase in rise and slight decrease in neighbouring 
backbone-backbone distance measured as salt concentra¬ 
tion is decreased (the duplex radius remains approxi¬ 
mately constant). Although there is not much experi¬ 
mental evidence to compare this with, there is some indi¬ 
cation that the pitch is roughly constant for low salt con¬ 


centrations (0.162 M and below).I^As mentioned earlier, 
the pitch is chosen (by modifying the bonded neighbour 
backbone-backbone interaction) so that the global twist 
of origami structures agrees with experimental measure¬ 
ments. Specifically, we set the backbone-backbone inter¬ 
action so that the helix bundle designed to have no global 
twist has no global twist in the model at [Na“'']=0.5 M. 
This results in a pitch of roughly 10.55 bp/turn at this 
salt concentration, compared to 10.36 bp/turn in the 
original model, and experimental val ues of around 10.45 
suggested by cyclisation experimentsalbeit in the 
presence of some divalent salt. 

The thermodynamics of duplex formation are shown in 
Fig.0 The transition width for the yield curve of a 10-bp 
duplex at 0.5 M [Na"*"] in oxDNA2 is largely unchanged 
from the original oxDNA, and the transition widths de¬ 
pend weakly if at all on salt in oxDNA2. The free-energy 
profiles for duplex formation in oxDNA2 show that the 
free-energy cost of forming the first base pair decreases 
with increasing salt, presumably due to the reduced ener¬ 
getic cost of bringing the two single strands close together 
as the electrostatic repulsion between backbone sites be¬ 
comes more short-ranged. The slope of the bound region 
of the free-energy profile also decreases with increasing 
salt, indicating a reduced free-energy cost of forming sub¬ 
sequent base pairs. 

In Fig. the duplex melting temperatures as a func¬ 
tion of salt and for different duplex lengths in oxDNA2 
are compared to the results from the SantaLucia model, 
to which oxDNA2 was parameterised. As expected, we 
find good agreement. 

The mechanical properties of the model double- 
stranded DNA as a function of salt, specifically the per¬ 
sistence length and torsional stiffness, are shown in Fig.|^ 
The persistence length was calculated by computing the 
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FIG. 7. The thermodynamics of formation of a 10-bp du¬ 
plex with oxDNA2. (a) Duplex yield for oxDNA2 at different 
salt concentrations, (b) Comparison of dnplex yield profiles 
at [Na+]=0.5M as predicted by oxDNA2, the original oxDNA 
model and the SantaLucia model.l^ (c) The free-energy profile 
for duplex formation at different salt concentrations, at the 
melting temperature of a 10-bp duplex at [Na"'"]=0.25 M. The 
duplex was simulated for at least 3 x 10^® virtual move Monte 
CarlcP^ (VMMC) steps at each salt concentration using um¬ 
brella sampling at a monomer concentration of 3.3 x 10“"^M. 


correlation of the helix axis as described in Ref. HZl The 
effective torsional stiffness, CeS, was computed from the 
linear regime of the torque response curve of a 60-bp du¬ 
plex under tension (see Appendix B2 for details). The 
simulations were carried out under a linear force of 30 pN, 
a high force regime where we expe ct th at CeS approxi¬ 
mates the true torsional stiffness Co The persistence 
length at [Na+]=0.5 M, about 123bp, is very slightly 
lower than for the original model, 125 bp, and the per¬ 
sistence length at different salt concentrations is consis¬ 
tent with t he rat her broad range of values suggested by 
experiment !II The decrease in persistence length with 
increasing salt is expected due to the decrease in repul¬ 
sion between the duplex’s backbone sites, which makes 
the duplex less stiff, in agreement with experiment The 
slight decrease in torsional stiffness with increasing salt 
can be rationalised in the same way. Regardless of the 


FIG. 8. The melting temperature as a function of salt con¬ 
centration for duplexes of different lengths simulated with 
oxDNA2 (crosses with solid lines) and as predicted by the San¬ 
taLucia model (crosses with dotted lines) for a DNA strand 
concentration of 3.3 x 10“^ M. Each duplex was simulated 
for roughly 4 x 10® VMMC steps at a salt concentration of 
0.25 M, and the melting temperatures at 0.15 M and 0.5 M 
were computed by single histogram reweighting.^ 


salt concentration, the torsional stiffness measured for 
oxDNA2 (Ceff « 380—400 fjfm) is lower than for the orig¬ 
inal modeP^ (Co « 473fJfm). Single-molecule twisting 
experiments give a value of Ceff « 410fJfm for a pulling 
force of 3.5 pN^and a pulling force of 45 pN,l^at around 
[Na“'"]=0.1 M. There is limited experimental data on the 
salt dependence of torsional stiffness in DNA; however 
the torsional stiffness has been reported to be roughly 
constant in the range [Na+]=0-0.162 M.S3 

The effects of various motifs on duplex melting tem¬ 
peratures in oxDNA2 are compared to results from the 
original oxDNA and the SantaLucia model in Table |TT] 
These effects are either barely changed in oxDNA2, or 
are now closer to the SantaLucia values than the original 
oxDNA values were. 


The melting temperatures of hairpins at 0.5 M [Na+] 
are shown in Fig. |10| The hairpin melting temperatures 
are lower than for the original oxDNA by about 2K, 
which in turn had hairpin melting temperatures about 
3 K lower than those predicted by SantaLucia. The gap 
between the oxDNA2 hairpin melting temperatures and 
those given by the SantaLucia model widens somewhat 
as the salt concentration is lowered (Fig. 11), indicating 
a stronger salt dependence in oxDNA2, with a typical 
underestimate of around 7K at a salt concentration of 
100 mM for a relatively short loop length of 6 bases. For 
longer loops, the difference in predicted melting temper¬ 
atures between oxDNA2 and SantaLucia widens further. 
This difference is unsurprising, as in oxDNA2 (and physi¬ 
cal DNA) ssDNA becomes stiffer at lower salt concentra¬ 
tions, making the formation of a hairpin less favorable,!^ 
whereas in the SantaLucia model the loops’ contribution 
to hairpin stability is salt- independent.!^ 
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FIG. 9. (a) The persistence length and (b) the torsional stiff¬ 
ness of duplex DNA in oxDNA2 as a function of salt con¬ 
centration. For the persistence length, a 60-bp duplex was 
simulated for at least 3 x 10® MD steps for each salt concen¬ 
tration. The method for computing the torsional stiffness is 
given in Appendix |B 2[ Error bars for the simulation results 
show the standard error on the mean given by averaging over 
either 10 independent estimates (for the persistence length) or 
5 independent estimates (for the torsional stiffness) for each 
data point. 



Number of bases in the loop 


FIG. 10. The melting temperature as a function of loop length 
for hairpins simulated with oxDNA2 (crosses with solid lines), 
the original oxDNA (crosses with dashed lines), and according 
to the SantaLucia model (crosses with dotted lines), at 0.5M 
[Na®"]. For the oxDNA2 results, the hairpins were simulated 
at each salt for at least 1.8 x 10®® VMMC steps. 


It seems plausible that oxDNA2’s performance is bet¬ 
ter than implied by the salt-independent loop contribu¬ 
tion in the SantaLucia model. In Fig. [T^ we consider 
hairpins with a short 5-bp stem and a long 31-base loop, 
making the hairpin thermodynamics particularly sensi¬ 
tive to the change in loop stiffness with salt. The hair¬ 
pin stability (i.e. the free energy difference between the 
bound and unbound state) as predicted by oxDNA2 and 
the SantaLucia model is compared to experimental data 



FIG. 11. The melting temperature as a function of salt 
concentration of hairpins with 6-base loops simulated with 
oxDNA2 (crosses with solid lines) and according to the San¬ 
taLucia model (crosses with dotted lines). For the oxDNA2 
results, the hairpins were simulated at each salt for at least 
1.5 X 10®° VMMG steps. 


determined using FRET (experimental details in Ap¬ 
pendix]^. As expected we find the hairpin stabilities for 
oxDNA2 have a much stronger salt dependence (steeper 
gradient) than those of SantaLucia (shallow gradient). 
At higher salt the oxDNA2 results show a similar slope 
to the experimental curve, but at lower salt (0.2 M and 
below) oxDNA2 shows a stronger destabilization of the 
hairpins with decreasing salt (a steeper gradient) than 
experiment. Thus, although oxDNA2 does a reasonable 
job of capturing these physical effects, the comparison 
suggests that the single strands may be experiencing too 
much repulsion at low salt. We also note that the San¬ 
taLucia model is unable to predict the large difference in 
stability of the two hairpins due to their different loop 
sequences, because the SantaLucia model is insensitive 
to the loop sequence except for the two bases adjacent to 
the stem. 


We examine the flexibility of single strands as a func¬ 
tion of salt in oxDNA2 by measuring the radius of gy¬ 
ration, Rg, of a single strand of DNA at different salt 
concentrations, sequences and strand lengths, and com¬ 
paring to experimental results (Fig. 13). The oxDNA2 
model reproduces the overall trend of increasing Rg with 
salt, and is in agreement with the somewhat noisy exper¬ 
imental data. We note that oxDNA2 is able to capture 
two important effects: the greater stiffness (and hence 
larger Rg) of the more strongly-stacking poly-A strands; 
and the greater salt dependence of the poly-T strands’ 
Rg compared to that of the poly-A strands, caused by 
the weaker stacking of the poly-T strands, which means 
that electrostatic repulsion makes a greater relative con¬ 
tribution to its stiffness. 
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_at^/k_ 

Motif oxDNA2 oxDNA SantaLucia 


Bulge (1-base bulge in 8-bp duplex) 

-18.7 

-17.98 

-23.19 

Bulge (2-base bulge in 8-bp duplex) 

-27.4 

-23.92 

-26.73 

Terminal mismatch (1-bp mismatch in 5-bp duplex) 

-b6.5 

-16.71 

-b8.6 

Internal mismatch (2-bp mismatch in 8-bp duplex) 

-15.9 

-15.77 

-13.99 


TABLE II. The effect of introducing various motifs on the duplex melting temperature at a monomer concentration of 3.3 x 
10“^ M. ATm is the difference between the melting temperature of the structure with the motif and a duplex consisting of the 
same number of complementary base pairs as the motif structure. For bulges and internal mismatches, the motif was placed 
at the centre of the duplex. The simulations of the bulges, terminal mismatches and internal mismatches were run for at 
least 6 X 10^°, 8.5 x 10^°, and 5.6 x 10^° VMMC steps respectively. All oxDNA2 results were obtained from simulations at a 
concentration of 4.2 x 10“® M which were extrapolated to a concentration of 3.3 x 10“^^ M. 
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FIG. 12. The free energy difference between the bound and 
unbound state for two hairpins as measured in experiment us¬ 
ing FRET (dashed lines), and as predicted by oxDNA2 (solid 
lines) and the SantaLucia model (dotted lines). Results for 
two different hairpins are presented: the A31-VW (red, circu¬ 
lar markers) hairpin and the T31-VW hairpin (blue, square 
ma rker s). The sequences for these hairpins are given in Ta¬ 
ble |A4] The oxDNA2 simulations were run for 10^° VMMC 
steps for each hairpin, with Qeff = 0 and a; = 0 (with x de¬ 
fined in Eq. |10| | , and the results (computed for the standard 
oxDNA2 values for geff and x) were calculated using ther¬ 
modynamic integration as described in Section |VI| and Ap¬ 
pendix |B 3| 


IX. CONCLUSION 

OxDNA has been applied to a wider variety of sys¬ 
tems than any other coarse-grained model of DNA. The 
modifications and extensions to oxDNA presented here 
open up a variety of new potential applications for 
the model. With the introduction of an explicit salt- 
dependent term in the potential, the model can be used to 
simulate systems under physiological conditions, and to 
investigate the salt-dependent behaviour of DNA. Also, 
the parametrization presented here allows a quantitative 
comparison to experiments run in a wide range of salt 
concentrations rather than just in the high-salt limit. 


FIG. 13. The radius of gyration {Rg) for ssDNA as a function 
of ssDNA length and for different salt concentrations, for (a) a 
poly-A ssDNA and (b) a poly-T ssDNA, from oxDNA2 (solid 
lines with circles) and from experiments due to Sim et 
(dashed lines with squares). The salt concentrations shown 
are 0.225M (blue lines), 0.525M (red lines), 1.025M (green 
lines). For the oxDNA2 results, the strands were simulated 
for between 2 x 10® and 5 x 10^® VMMG steps at each strand 
length and salt concentration. 


The introduction of major-minor grooving adds detail to 
the model, and, combined with other small modifications, 
allows the use of oxDNA2 to accurately characterise the 
structural properties of large DNA nanostructures, such 
as DNA origami. In particular the helical pitch and the 
twist angles at nicks and junctions were fine-tuned to 
obtain a correct global twist for a test-case 3D origami 
structure. In the absence of definitive experimental val¬ 
ues for these individual parameters we chose a combi¬ 
nation that we deemed physically reasonable and that 
produces equilibrium structures that compare well with 
experimental ones. Finally, the sequence dependence in 
the model has been extended by introducing different in¬ 
teraction strengths for the AA and TT stacking. This 
change will be particularly useful for studying the ef¬ 
fects of stacking in single strands or single-stranded sec¬ 
tions (e.g. hairpin loops) as poly-A and poly-T sections 
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are often used as paradigmatic examples of strongly and 
weakly stacking sequences, respectively, and for mod¬ 
elling DNA nanostructures where poly-T loops are often 
used as flexible linkers. 

There are some limitations associated with the new 
features of the model. The Debye-Hiickel treatment of 
the electrostatics is perhaps an oversimplification, and 
we should not expect it to capture all of the complexi¬ 
ties associated with the electrostatics for DNA, but it is 
not straightforward to think of a different approach that 
would be consistent with our level of coarse graining. In 
particular, the fact that hairpin stability is reduced faster 
than in experiment suggests that single strands may ex¬ 
perience too much repulsion in oxDNA2 at low salt, and 
care should be taken in making predictions based on our 
models at low salt concentration if the system under in¬ 
vestigation depends crucially on the thermodynamics of 
long single-stranded sections. 

We are currently exploiting the improvements in the 
model’s structural prediction to carry out investigations 
on DNA origami structures as well as on structures com¬ 
posed of multi-arm tile^^. At the same time, we are 
studying the salt dependent thermodynamics of a diverse 
set of systems, which would not have been possible with¬ 
out the introduction of salt dependence into the poten¬ 
tial. The introduction of different strengths for the AA 
and TT stacking in the model allows us to better cap¬ 
ture these effects, and in addition gives us more accurate 
sequence-dependent hairpin thermodynamics and kinet¬ 
ics, which are currently being exploited to further study 
hairpins. We note that a simulation code implementing 
both the original oxDNA model and the new oxDNA2 
model, including an implementation using GPUs,l^ is 
available for download from dna. physics .ox.ac.uk. 
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Appendix A: Modifications to the oxDNA potential in oxDNA2 

1. Introducing an electrostatic term into the model potential 

We introduce an explicit electrostatic term into the model using Debye-Hiickel theory. To aid computational 
efficiency when simulating the model, we ensure that the interaction goes to zero within a finite distance Ccut.DH, 
and we introduce a quadratic smoothing function so that the interaction goes to zero smoothly. The form of the 
electrostatic term is then 


Kslectrostatic (?■ jT,!) 


V;^.ooth(r‘’-^T, I) 
0 


if ^smooth,DH ^ ^ : 

if ^cut,DH ^ r ^ r'smooth.DH; 

otherwise, 


(Al) 


where T is the temperature, I is the salt concentration, and Tsmooth.DH, the distance above which smoothing is 
introduced, is chosen to be equal to 3Adh (Adh is defined below). The Debye-Hiickel-like term, T’, d), is 

given by 


vbH(A^^,J) 


(geffc)^ exp Adh(T, I)} 

dTreoCr 


(A2) 


where geff is the effective charge (which is equal to 1 in Debye-Hiickel theory but which we choose to be 0.815 after 
fitting to experimental data), e is the charge of an electron, cq is the permittivity of a vacuum, is the relative 
permittivity of water (which we choose to be 80), and is the distance between the pair of interacting backbone 
sites. Adh is the Debye screening length, and is given by 


•^dh(T", I) 


eoerkBT 
2NAe^I ’ 


(A3) 


where fcs is Boltzmann’s constant and Na is Avogadro’s number. The quadratic smoothing term, K 5 niooth(c'’’ 
is given by 


k^smooth(r'’■^ T, I) = (A4) 

and B and rcut.DH are constants chosen such that Veiectrostatic is smooth and differentiable. Note that the above 
definition implies that the cutoff radius for the Debye-Hiickel term, rcut.DH, depends on the salt concentration I and 
the temperature T. 


2 . Modifying the average pitch of the model duplex 

The model duplex pitch can be controlled by modifying the position of the minimum in the FENE bonded back¬ 
bone potential, f^rj^g^ckbonej which has a small effect on the thermodynamics (which was refitted using the histogram 
reweighting method described in Section H C 2 of Ref. ICT) and little effect on the other properties of the model DNA. 
Eor oxDNA2, was set to give a global twist of zero for the N-type helix bundle due to Dietz et a/.j^^the new 


and old values for are given in Table |A1 


3. Modifying the average coaxial stacking term in the oxDNA potential 


We modified the coaxial stacking term in the oxDNA potential in order to change the twist angle across a nick 
and across an origami junction, while changing the thermodynamic properties of a nicked duplex and other motifs as 
littl e as possible. This was achieved by moving the position of the minimum with respect to the angle 9i (shown in 
Fig. Al) of the coaxial stacking term, and by increasing the well depth, fccoax, to make the minimum more 

energetically favourable. Increasing the well depth was necessary because it is less favourable for a pair of coaxially 
stacking nucleotides in a nicked duplex to stack with this tighter angle; a stronger overall coaxial stacking interaction 
compensates for this effect. The changes to the model are shown in Table [AT] 

The altered optimal angle i requires an alteration of the modulation function involving 9i for coaxial stacking. 
In the original oxDNA, the 9i modulation is given by / 4 (di) -|- — 9i), where /4 is defined in Ref. HSl The second 
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FIG. Al. Schematic showing the angle 9i between two nucleotides. Each nucleotide is represented by three circles joined by 
a line; the large solid circles represent the backbone sites, while the small solid circles represent the stacking (closer to the 
backbone) and hydrogen-bonding (at the end of the nucleotide) sites. In this paper we modify the term dependent on 6\ for 
the coaxial stacking interaction in the oxDNA model, in order to reduce the overtwist across nicks and junctions. 


Parameter 

oxDNA 

oxDNA2 

7/degrees 

0 

20 

Qefi 

0 

0.815 

Cr 

- 

80 

^coax 

46 

58.5 

^0 

^coaXjl 

TT — 0.60 

TT - 0.25 

-^coax. 1 

- 

40 

-^coax, 1 

- 

TT - 0.025 

ehb 

1.077 

1.0678 

€stack 

1.3448 -f 2.6568fcBr 1.3523 -f 2.6717/03^ 

X.f.0 

' backbone 

0.7525 

0.7564 


TABLE Al. Comparison of model parameters between oxDNA and oxDNA2. All other parameters are unchanged in oxDNA2. 
All angles except 7 are given in radians, all lengths are defined with respect to a reduced length scale (1 unit = 8.518A) and 
all energies are defined with respect to a reduced temperature {k^T = 0.1 corresponding to 300K). 7 is the angle between the 
line from the helix axis to the backbone site and the line from the helix axis to the stacking site (see Eig. [^, which defines 
the backbone site position in oxDNA2. qeff is the effective charge used for the Debye-Hiickel treatment introduced in the new 
model for the electrostatic interactions due to the negatively-charged phosphate groups in the DNA backbone, tr is the relative 
permittivity of water which is used for the Debye-Hiickel treatment, fccoax determines the overall strength of the coaxial stacking 
interaction while ^coax.i specifies the position of the minimum with respect to the angle 9i of the coaxial stacking potential. 
The values for Acoax.i and Hcoax.i are chosen to ensure that the coaxial stacking potential is smooth and differentiable, chb and 
Estack determine the well depths for the hydrogen bonding and stacking interactions respectively. gives the position 

of the minimum of the FENE bonded backbone potential. 


term was introduced to ensure that the potential is differentiable at 9i = tt, where the gradients of the two terms 
cancel. In oxDNA, this additional term had a very small effect on the numerical value of the potential for 61 < tt. 
Using the same trick in oxDNA2, however, is problematic, because ^coax.i is much closer to tt. The second term, 
/4(27r — 0i), would then be substantial for 0i < tt and would distort the modulating function in an undesired fashion. 
To avoid this, we replace the second term, Ja^^tt — 0i), with a new function fei0i)- 

^{0-Bf. ii9>B, 

0 otherwise. 

The f 4 { 0 i) + fii^TT — 0i) factor in the coaxial stacking term is then replaced by f 4 { 0 i) + /6(^i)- The values for Acoax,i 
and i?coax.i are chosen to ensure that the coaxial stacking potential is smooth and differentiable; they are specified in 

Table EH 

For oxDNA2, we made an additional change to the coaxial stacking potential: we removed the 
/ 5 (cos((/) 3 ),acoax. 3 ',cos((/) 3 )*oa,,) and /5(cos((/)4),acoax.4',cos((/)4)coax) terms. These terms had allowed only right- 
handed blunt-ended stacking, disallowing left-handed blunt-ended stacking; however, since the development of the 
original oxDNA model, this constraint has been deemed unnecessary, as the blunt-ended coaxial stacking interaction 
was found to be likely to be too weak compared to experiment,!^ and there is no experimental evidence indicating 
that left-handed blunt-ended coaxial stacking is not possible. 
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Measurement 

oxDNA 

oxDNA2 

SantaLucia 

AG‘nick_stack 

-4.3fcBT 

-4.4fcBT 

- 

AXrn, hairpin _stem,6mer 

-b9.9K 

-bl0.6K 

-bll.7K 

ATjn^hairpin-stemjTmer 

-b7.7K 

-b8.3K 

-t8.5K 

AXrn,hairpin_stem,8mer 

-b6.3K 

-b6.9K 

4 -6.4K 

AG" blunt ended stacking 

-b4.9feBr -b5.3fcBT 

- 


TABLE A2. Comparison of some thermodynamic quantities that depend on the coaxial stacking interaction in oxDNA and 
oxDNA2. AGnick_stack = AGstacked “ AGunstacked is the free-energy change upon stacking across a nick for a 20-bp duplex 
with a nick at the centre at 37° C. Arin,hairpin_stem,6mer = T'm,hairpin_stem,6mer “ Tm.emer gives the difference in melting tern - 
perature at a concentration of 3.3 x 10““^ M between a 6-bp duplex adjacent to a hairpin and a 6-bp duplex (see Fig. A2l, 
and ATin,hairpin_stem,7mer and Arm.hairpin_stem,8mer have analogous definitions. AGbiunt_ended_atacking is the fiee-energy change 
upon dimerisation of two 6-bp duplexes through blunt-ended coaxial stacking at a monomer concentration of 5.37 mM and at 
19.85° C. AGnick_atack was Computed from a simulation of 1.9 x 10^^ VMMC steps, each ATni,hairpin_atem,Xmer from simulations 
of 8 X 10^° VMMC steps, and AGbiunt_ended_atacking from a simulation of 2.4 x 10^^ VMMC steps. 



+ 


FIG. A2. A schematic showing the transition important for the calculation of Arin,hairpin_stem,Xmer (see Table [A2| , which is 

given by ATnijhairpin_steni,Xmer '— Tni,hairpin_stem,Xmer -7iii,xmer, rvhere 7in,hairpiii_stem,Xnier is the melting temperature of a 1C bp 

duplex adjacent to a hairpin (the schematic shows this system for the X = 4 case) and Tm.Xmer is the melting temperature of 
a X-bp duplex. 


Incorporating the changes described above, the final form for the oxDNA2 coaxial stacking term is 




coax_stack — 


/2((5ratack, fccoax, SrlZ^, (5r^oix)/4(6'4, Ocoax.d, Cax,4: 

^ (/4(^lj Ocoax.t) ^coax,li ^^coax,r) d” : -^coax, 15 .^coax,r ) ) 

^ (^*4(^5 ; ^coax,5; ^coax,5 5 ^^coax.s) ; ^coax,5 5 ^coax,5 ’ ^^coax.s)) 

^ (^* 4 (^ 6 ; acoax,6; ^coax,65 ^^coax.s) d~ /4(^ ^6; acoax,65 ^coax,6; ^^coax.s))* 


(A6) 


Table |A2| gives a comparison of the thermodynamic properties of the motifs relevant to the coaxial stacking term 
in oxDNA and oxDNA2. The free-energy change upon stacking across a nick, AGnick_stackj is slightly more negative 
in oxDNA2, indicating that a nicked duplex stacks slightly more strongly across the nick in oxDNA2 than in the 
original oxDNA. In both cases the stacking is stronger than is seen in experiment (—2.62 A:bT in Ref. IMl - a discussion 
comparing the value of AGnick_stack measured for the original oxDNA with experimental results is given in Ref. I28|l . 
The stabilising effect of a hairpin stem on duplex hybridisation (see Fig. |A2[ ATm hairpin_stem,Xmer) is found to be 
slightly stronger than for the original oxDNA, and is on average in better agreement with the stability predicted by 
the SantaLucia model. Finally, dimerisation of two duplexes through blunt-ended coaxial stacking is found to be less 
favourable in oxDNA2 than in oxDNA. This is despite the mitigating effect of removing the /s terms in the coaxial 
stacking interaction, which strengthens blunt-ended stacking. A previous studjEH has shown that oxDNA already 
underestimated the stability of blunt-ended coaxial stacking. 
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4. Modifying the hydrogen bonding and stacking interaction strengths 

Model duplexes were found to fray slightly less often when different widths for the major and minor duplex grooves 
were introduced compared to the equal-groove-width case, which led to a slight raising of the duplex melting temper¬ 
atures. In addition, the change to the position of the minimum in the FENE bonded backbone potential, ^?'backboneJ 
was found to destabilise the duplex state, slightly lowering the duplex melting temperatures. In the process of the 
reparameterisation of oxDNA to incorporate these two changes, the duplex melting temperatures were reset on the 
introduction of each change by making small modihcations to the stacking and hydrogen bonding strengths using the 
histogram reweighting method described in Section II C 2 of Ref. I6IL giving an agreement with experimental melting 
temperatures that was just as good as before. The new parameters are specified, and compared to the old parameters, 
in Table |A1[ while Fig. HI gives a comparison of the duplex thermodynamics with the old and new versions of the 
model. 


Appendix B: Simulation methods 

1. Measuring the global twist of an origami 

We measure the global twist of the origami due to Dietz et a/.^ in the following way: The origami (for example, 
the L-type, although this applies equally well to all three types) can be thought of as a set of antiparallel double 
helices on a hexagonal lattice joined together by crossovers. The overall shape of this structure is roughly that of a 
(possibly twisted) rectangular prism. To measure the twist, we consider the two faces of the rectangular prism which 
are opposite each other and approximately normal to the axes of the double helices. For each face, we define a vector 
along its long edge, and we superimpose these two vectors into a plane perpendicular to the average of the helix axes. 
We then define the global twist as the angle between these two vectors. 

To avoid end effects, the top and bottom rows of helices (3 helices in the top row and 3 helices in the bottom row) 
were excluded for this analysis, as were approximately 22 base pairs from each end of every double helix. 


2. Measuring the torsional stiffness of a duplex in the oxDNA model 


The torsional stiffness of a duplex in oxDNA2 was measured according to the scheme outlined by Matek et alW^ In 
summary, this method involves pulling a duplex while also twisting it, and measuring the torque required to impose 
a given twist. In this case, a 60-bp duplex was pulled at a constant force of 30 pN and held with virtual traps which 
imposed a helical pitch different from the equilibrium pitch. For each salt concentration, simulations were run with 
imposed pitches of 9.8, 10.1, 10.4, 10.7, 11.0, and 11.3 bp/turn, with each simulation run for at least 2.6 x 10® MD 
steps. The torque observed was plotted as a function of imposed twist (this is called a torque response curve), and 
the gradient of the linear region used to approximate the effective torsional stiffness. 

The torsional stiffness Ceft is computed as 


Ceff = 


AF Co 
Act 9n ’ 


(Bl) 


where AF is the change in torque, and Act the change in superhelical twist density, measured in the linear regime of 
the torque response curve, oq is the double helical rise for a relaxed duplex and 6 q is the twist angle across a base-pair 
step for a relaxed duplex, measured in radians. According to Moroz and Nelson,!^ the effective torsional stiffness Ces 
measured by twisting a duplex under tension can be related to the true torsional stiffness Cq by 


= 1 -^ 


kBT ] 

BoFj 


(B2) 


where F is the linear force and Bq is the bending stiffness. Our simulations are in the high-force regime, where CeS 
should be a good approximation for Cq. 


3. Details of the simulation protocols 

The simulation results presented in this paper were obtained using either a virtual move Monte Carlo (VMMC) 
algorithm or a molecular dynamics (MD) algorithm with a Brownian thermostat. Note that, where VMMC steps are 
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Parameter 




Simulation description 




Standard 

Bulge 

Nick stack. Hairpin stab. 

Blunt coax 

Hairpin SF C7torsionai 

Ip 

Helix bundle twist 

Algorithm 

VMMC 

VMMC 

VMMC 

VMMC 

VMMC 

VMMC 

MD 

MD 

MD 

^trans 

0.11 

0.1 

0.03 

0.15 

0.22 

0.11 

- 

- 

- 

^rot 

0.22 

0.2 

0.15 

0.15 

0.22 

0.22 

- 

- 

- 

max. cluster size 

- 

- 

- 

- 

12 

- 

- 

- 

- 

Verlet skin 

1 

1 

1 

1 

1 

1 

0.05 

0.1 

0.05 

^^step 

- 

- 

- 

- 

- 

- 

0.005 

0.005 

0.005 

diffusion coefficient 

- 

- 

- 

- 

- 

- 

2.5 

2.5 

2.5 

T 


Tm 

310.15K 

Tm 

293K 

295.6K 

296.15K 

296.15K 

296.15K 


TABLE A3. The results in this paper were computed using simulations of oxDNA2 with these parameters. The “standard” 
parameters were used for the following calculations: ATm for a terminal mismatch and ATm for an internal mismatch (Table|Il|; 
all hairpin Tm calculations (Fig. 10 and Fig.jHi; and all duplex Tm calculations (Fig.[^. The “bulge” parameters were used for 
the ATm calculations for a bulge motif (Table the “nicked stacking” parameters were used for the AGnick_stack calculation 

(Table [A^ , the “hairpin stem stab.” parameters were used for the ATm,hairpin_stem,Xmer calculations (Table A2|, the “blunt 
coax” parameters were used for the calculation of AGbiunt_ended_stacking (Table A2|, the “hairpin df”’ parameters were used 
for the initial free energy calculation for the thermodynamic integration to find the relative stability of the hairpins shown in 
Fig. 12 and to find the radius of gyration {Rg) of ssDNA (Fig. 131, the “Gtorsionai” parameters were used for the calculation 
of the torsional stiffness (Fig.j^b)), the “Ip" parameters were used for the calculation of the persistence length (Fig. [^a)) and 
the pitch and rise (Fig. [^, and the “helix bundle twist” parameters were used for the calculation of the global twist of the 
helix bundle origami structures (Table |^. dtrans and 5rot are the standard deviation of the normally distributed move size for 
the translational and rotational VMMC moves respectively, max. cluster size is the size of the biggest cluster that VMMC is 
permitted to choose for a trial move, Verlet skin is the thickness of the Verlet skin used for Verlet lists, Ststep is the size of 
the simulation time step in MD, diffusion coefficient gives the diffusion coefficient used for a single nucleotide, and T is the 
temperature of the heat bath {Tm indicates that simulations were run at or close to the melting temperature for that system). 
All numbers are given in oxDNA units (1 oxDNA length unit = 0.8518 nm, 1 oxDNA time unit = 3.03 ps), unless given explicitly 
(angles are quoted in radians). 


reported in the main text, they refer to the number of accepted moves during the simulation. Details of the simulation 
parameters used are given in Table |A3[ 

For the duplex free-energy profiles, yields and melting temperatures (Fig. and Fig. we used umbrella sampling 
with an order parameter defined in terms of the number of native base pairs (0 to n where n is the number of base 
pairs in a fully formed duplex). A base pair was defined as being formed if the potential’s hydrogen-bon ding energy 


10 


Fig. 


11 


term for that pair was lower than —0.596 kcalmol ^. For the hairpin melting temperature calculations (Fig. 
and Fig. 12) we used an analogous order parameter for the native base pairs in the hairpin stem. 

The oxDNA2 results reported in Fig.were computed using thermodynamic integration as described in Section [Vl| 
For each oxDNA2 data point shown in the figure, first an umbrella sampling simulation of each hairpin (sequences 
shown in Table A4) was carried out for a version of th e model with x = Q, where x is given by Eq. 10 and Qeff = 0 in 
order to find the reference free energy in Eq. 11). The simulations were run for 6 x 10^^ VMMC steps and used 


the parameters specified in Table |A3| The integrals in Eq. IT were discretised into 15 gridpoints for the integral in 
ggff and 10 gridpoints for the integral in x. Then simulations to compute these integrals were run for roughly 6 x 10^ 
VMMC steps for each grid point, so that each data point in the figure required 25 grid point simulations, with one 
umbrella sampling simulation per hairpin. 


Appendix C: Experimental Methods 
1. DNA Hairpin Design 

The hairpin constructs consist of a loop with different lengths (21 or 31 nucleotides) and sequences (adenosine or 


thymine), and six base pairs in the stem, as shown in Table A4 


A 35-bp long, double-stranded DNA was used to keep the hairpin section away from the coverslip and the origami 
surfaces. The hairpins were labeled with donor and acceptor fluorophores (ATTO-550 and ATTO-647N, respectively), 
positioned such that the hairpin open state yielded low FRET values {E) and the closed state yielded high E values 
(Figure]^. 
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Name Hairpin Sequence 

T 21 -W 5 '-TGGGTT-(T) 2 i-AACCCA- 3 ' 
T 31 -VW 5 '-TGGATT-(T) 3 i-AATCCA- 3 ' 
A 21 -W 5 '-TGGGTT-(A) 2 i-AACGGA- 3 ' 
A 31 -VW 5 '-TGGATT-(A) 3 i-AATGCA- 3 ' 


TABLE A4. DNA hairpins used for oxDNA validation. Note that the experimental system included an additional dangling 
duplex attached to the hairpin (denoted “dsDNA part” in Fig. A31. For the oxDNA2 simulations, we used a truncated (13-bp) 
duplex in order to capture any effects due to the duplex interfering with the hairpin while avoiding the larger computational 
cost of simulating the full duplex. 


Loop 


Stem 



dsDNA part 


XXXXXA—TT--TTCTTCACAAACCAGTCCAAACTATCACAAACTTA-3' 
YYYYYB-k-sl-AAGAAGTGTfTlTGGTCAGGTTTGATAGTGTTTGAAT-S' 


r^O-55^ 



FIG. A3. Design of the DNA hairpin constructs used to parameterise the oxDNA model. The hairpin sequences and labeling 
positions were designed to minimize possible interfering interactions between the fluorophores when the hairpins were closed, 
and the TT spacer was introduced to minimize interactions between the stem and the duplex. 


2. Annealing Procedures and Experimental Setup 

The top and bottom strands were annealed at 94°C (1.5 p-M in 10 pL TE-NaCl buffer [10 mM Tris, pH 8.0, ImM 
EDTA], and 100 mM NaCl) and then gradually cooled (30 min) to room temperature (using PCR). 

The measurements were performed on a diffusion-based single-molecule FRET-ALEX setup, as described 
elsewhere.l^ In brief, a green CW laser beam (532 nm, CL532-025-L, Crystal Laser, Reno, NV, USA) was 
aligned/misaligned into a single-mode fiber using an acousto-optic modulator (AOM; R23080-2-LTD, Neos Tech¬ 
nologies, Melbourne, EL, USA), alternating with a red diode laser (640nm, 1069417, Coherent, Auburn, CA, USA) 
that was electronically switched on/off. The AOM and the red laser were computer controlled with a 12.4-ps on-time, 
a 12.6-ps off-time, and a phase shift of 12.5 ps. The intensity rise and fall times were less than 50 ns, and there was 
no time overlap between the lasers. The laser beams were combined by a dichroic mirror (Z532RDC, Chroma, Bel¬ 
lows Falls, VT, USA) and coupled into a single-mode fiber (P1-460A-FC-2, Thorlabs, Newton, NJ, USA). The laser 
intensities were tuned such that the doubly labeled species would yield a stoichiometry, S, of roughly 0.5 (90 pW for 
the green laser and 35 pW for the red laser, measured after the fiber while alternating). After collimation (objective 
PLCNlO/0.25, Olympus America, Melville, NY, USA), the combined green and red beams were introduced into a 
commercial inverted microscope (1X71, Olympus America) and focused about 70 pm inside the sample solution with 
a water-immersion objective (NA 1.2, 60x, Olympus America). The emitted fluorescence was separated from the 
excitation light by a dichroic mirror (ZT532/638RPC, Chroma), focused into a 100-pm pinhole (PIOOS, Thorlabs), 
re-collimated, split by a second dichroic mirror (FF650-Di01, Semrock, Lake Forest, IL, USA), passed through fil¬ 
ters (band-pass filter, Semrock FFOl-580/60, for the donor channel and a long-pass filter, Semrock BLP01-635R, for 
the acceptor channel), and focused into two single-photon avalanche photodiodes (SPAD; SPCM-AQRH-13, Perkin- 
Elmer Optoelectronics, Fremont, CA, USA). The TTL signals of the two SPADs were recorded as a function of time 
by a 12.5-ns resolution counting board (PCI-6602, National Instruments, Austin, TX, USA) and in-house prepared 
LabView acquisition software. 


3. Measurement Solution and Conditions 

Measurements were carried out on diluted samples (40 pL of 1 pM hairpins) that were placed on a coverslip and 
sealed with silicone and an upper coverslip. The bottom coverslip was KOH-treated to prevent sticking to the surface, 
sonicated for 15 min in 1 M KOH solution, thoroughly washed with distilled water, and dried in air. The measurement 

























21 



FIG. A4. Typical data for hairpin two-dimensional i5/5'-histogram and E- and S- one-dimensional histograms. The specific 
data refers to A31-VW measured in buffer containing 100 mM NaCl. Donor-only (green) and acceptor-only (red) populations 
are rejected by S-ratio. The populations with S « 0.5 (inside black rectangle) were used to obtain the hnal U-histogram.^^ 


buffer comprised 10 mM Tris pH 8.0, ImM EDTA, 10 pg/mL bovine serum albumin (BSA, Sigma-Aldrich) to reduce 
sample sticking, 2 mM Trolox (Sigma-Aldrich) to reduce photo-bleaching and photo-blinking of the dyes, and indicated 
NaCl concentrations. Data were collected for over 30 minutes, and temperature was maintained at 22.5° C. 


4. Data analysis 

Data analysis was performed with the in-house written LabView software as described elsewhere.!^ The beginnings 
and endings of bursts were determined by the all-photons-burst-search (APBS, parameters: L = 200, M = 10, and 
T = 500 ps. For each burst, E and S were calculated according, respectively, binned (0.01 bin size), and plotted 
on one-dimensional E and S histograms and on a two-dimensional E/S histogram.!^ Donor-only and acceptor-only 
populations were rejected by ALEX . Th e Obtained if-histogram showed two distinct peaks and was assumed to have 
quasi two-state dynamic^ (Figure [A4| . 

The opening and c losing rates were extracted from the shape of the E-histogram using dynamic Probability Distri¬ 
bution Analysis 


5. Experimental Data 

The results of the hairpin kinetics experiments are shown in Fig. |A5| 
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FIG. A5. The free energy difference between the bound and unbound state {5F) experimentally determined for the four 
hairpins. 




